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Abstract 

We consider the use of language models 
whose size and accuracy are intermedi- 
ate between different order n-gram models. 
Two types of models are studied in partic- 
ular. Aggregate Markov models are class- 
based bigram models in which the map- 
ping from words to classes is probabilis- 
tic. Mixed-order Markov models combine 
bigram models whose predictions are con- 
ditioned on different words. Both types 
of models are trained by Expectation- 
Maximization (EM) algorithms for maxi- 
mum likelihood estimation. We examine 
smoothing procedures in which these mod- 
els are interposed between different order 
n-grams. This is found to significantly re- 
duce the perplexity of unseen word combi- 
nations. 

1 Introduction 

The purpose of a statistical language model is to as- 
sign high probabilities to likely word sequences and 
low probabilities to unlikely ones. The challenge 
here arises from the combinatorially large number 
of possibilities, only a fraction of which can ever be 
observed. In general, language models must learn 
to recognize word sequences that are functionally 
similar but lexically distinct. The learning problem, 
one of generalizing from sparse data, is particularly 
acute for large-sized vocabularies (Jelinek, Mercer, 
and Roukos, 1992). 

The simplest models of natural language are n- 
gram Markov models. In these models, the prob- 
ability of each word depends on the n — 1 words 
that precede it. The problems in estimating ro- 
bust models of this form are well-documented. The 
number of parameters — or transition probabilities — 
scales as V n , where V is the vocabulary size. For 



typical models (e.g., n = 3, V — 10 4 ), this num- 
ber exceeds by many orders of magnitude the total 
number of words in any feasible training corpus. 

The transition probabilities in n-gram models arc 
estimated from the counts of word combinations in 
the training corpus. Maximum likelihood (ML) esti- 
mation leads to zero- valued probabilities for unseen 
n-grams. In practice, one adjusts or smoothes (Chen 
and Goodman, 1996) the ML estimates so that 
the language model can generalize to new phrases. 
Smoothing can be done in many ways — for example, 
by introducing artificial counts, backing off to lower- 
order models (Katz, 1987), or combining models by 
interpolation (Jelinek and Mercer, 1980). 

Often a great deal of information is lost in the 
smoothing procedure. This is due to the great dis- 
crepancy between n-gram models of different order. 
The goal of this paper is to investigate models that 
are intermediate, in both size and accuracy, between 
different order n-gram models. We show that such 
models can "intervene" between different order n- 
grams in the smoothing procedure. Experimentally, 
we find that this significantly reduces the perplexity 
of unseen word combinations. 

The language models in this paper were evalu- 
ated on the ARPA North American Business News 
(NAB) corpus. All our experiments used a vo- 
cabulary of sixty-thousand words, including tokens 
for punctuation, sentence boundaries, and an un- 
known word token standing for all out-of- vocabulary 
words. The training data consisted of approxi- 
mately 78 million words (three million sentences); 
the test data, 13 million words (one-half million 
sentences). All sentences were drawn randomly 
without replacement from the NAB corpus. All 
perplexity figures given in the paper are com- 
puted by combining sentence probabilities; the prob- 
ability of sentence wqWi ■ ■ ■ w n w n +i is given by 
P(wi\w Q ■ ■■ Wi-i), where w and w n+1 are 
the start- and end-of-sentence markers, respectively. 



Though not reported below, we also confirmed that 
the results did not vary significantly for different ran- 
domly drawn test sets of the same size. 

The organization of this paper is as follows. 
In Section [^, we examine aggregate Markov mod- 
els, or class-based bigram models (Brown et al., 
1992) in which the mapping from words to classes 
is probabilistic. We describe an iterative algo- 
rithm for discovering "soft" word classes, based on 
the Expectation-Maximization (EM) procedure for 
maximum likelihood estimation (Dempster, Laird, 
and Rubin, 1977). Several features make this algo- 
rithm attractive for large- vocabulary language mod- 
eling: it has no tuning parameters, converges mono- 
tonically in the log-likelihood, and handles proba- 
bilistic constraints in a natural way. The number 
of classes, C, can be small or large depending on 
the constraints of the modeler. Varying the number 
of classes leads to models that are intermediate be- 
tween unigram (C = 1) and bigram (C = V) models. 

In Section ||, we examine another sort of "inter- 
mediate" model, one that arises from combinations 
of non-adjacent words. Language models using such 
combinations have been proposed by Huang et al. 
(1993), Ney, Essen, and Kncser (1994), and Rosen- 
feld (1996), among others. We consider specifically 
the skip-k transition matrices, M(w t -k,Wt)i whose 
predictions are conditioned on the kth previous word 
in the sentence. (The value of k determines how 
many words one "skips" back to make the predic- 
tion.) These predictions, conditioned on only a sin- 
gle previous word in the sentence, are inherently 
weaker than those conditioned on all k previous 
words. Nevertheless, by combining several predic- 
tions of this form (for different values of k), we can 
create a model that is intermediate in size and ac- 
curacy between bigram and trigram models. 

Mixed-order Markov models express the predic- 
tions P{wt\wt-\,wt-2, ■ ■ ■ ,Wt-m) as a convex com- 
bination of skip-A: transition matrices, M(wt-k,Wt)- 
We derive an EM algorithm to learn the mixing co- 
efficients, as well as the elements of the transition 
matrices. The number of transition probabilities in 
these models scales as mV 2 , as opposed to v m+1 . 
Mixed-order models are not as powerful as trigram 
models, but they can make much stronger predic- 
tions than bigram models. The reason is that quite 
often the immediately preceding word has less pre- 
dictive value than earlier words in the same sentence. 

In Section |], we use aggregate and mixed-order 
models to improve the probability estimates from 
n-grams. This is done by interposing these models 
between different order n-grams in the smoothing 
procedure. We compare our results to a baseline tri- 



gram model that backs off to bigram and unigram 
models. The use of "intermediate" models is found 
to reduce the perplexity of unseen word combina- 
tions by over 50%. 

In Section |[ we discuss some extensions to these 
models and some open problems for future research. 
We conclude that aggregate and mixed-order models 
provide a compelling alternative to language models 
based exclusively on n-grams. 

2 Aggregate Markov models 

In this section we consider how to construct class- 
based bigram models (Brown et al., 1992). The 
problem is naturally formulated as one of hidden 
variable density estimation. Let P(c\wi) denote the 
probability that word w\ is mapped into class c. 
Likewise, let P(w2\c) denote the probability that 
words in class c are followed by the word W2- The 
class-based bigram model predicts that word u>i is 
followed by word u>2 with probability 

c 

P(W 2 \W 1 ) = J2 P ( W 2\C)P(C\W 1 ), (1) 
c=l 

where C is the total number of classes. The hidden 
variable in this problem is the class label c, which 
is unknown for each word w\. Note that eq. ([j]) 
represents the V 2 elements of the transition matrix 
P{w2\wi) in terms of the 2CV elements of P(w2\c) 
and P{c\wi). 

The Expectation-Maximization (EM) algorithm 
(Dempster, Laird, and Rubin, 1977) is an iterative 
procedure for estimating the parameters of hidden 
variable models. Each iteration consists of two steps: 
an E-step which computes statistics over the hidden 
variables, and an M-step which updates the param- 
eters to reflect these statistics. 

The EM algorithm for aggregate Markov models 
is particularly simple. The E-step is to compute, for 
each bigram w\W2 in the training set, the posterior 
probability 



P{c\w\, w 2 ) 



P(w2\c)P(c\wi) 

j2 cr p(w 2 W)P(c>\ Wl ) 



(2) 



Eq. (^j gives the probability that word W\ was as- 
signed to class c, based on the observation that it 
was followed by word W2- The M-step uses these 
posterior probabilities to re-estimate the model pa- 
rameters. The updates for aggregate Markov models 
are: 

T, w N(w 1 ,w)P(c\w 1 ,w) 



P(c\ Wl ) 
P{w 2 \c) 



Y. WC ' N(wi,w)P(cf\wi,w) ' 

J2w N ( w ' w ^) p ( c \ w ' w 2) 

Xw N ( w > w')P(c\w, w') ' 



(3) 
(4) 
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test 


1 


964.7 


964.9 


2 


771.2 


772.2 


4 


541.9 


543.6 


8 


399.5 


401.5 


16 


328.8 


331.8 


32 


278.9 


283.2 


V 


123.6 





Table 1: Perplexities of aggregate Markov models on 
the training and test sets; C is the number of classes. 
The case C = 1 corresponds to a ML unigram model; 
C = V, to a ML bigram model. 








winning assignment probability 



Figure 2: Histogram of the winning assignment 
probabilities, max c P(c|w), for the three hundred 
most commonly occurring words. 

where N(wi,W2) denotes the number of counts of 
W1W2 in the training set. These updates are guar- 
anteed to increase the overall log-likelihood, 

1= N(w 1 ,w 2 )lnP(w 2 \w 1 ), (5) 
W1W2 

at each iteration. In general, they converge to a local 
(though not global) maximum of the log-likelihood. 
The perplexity V* is related to the log-likelihood by 
V* = e~ e / N , where N is the total number of words 
processed. 

Though several algorithms (Brown et al., 1992; 
Pereira, Tishby, and Lee, 1993) have been proposed 
for performing the decomposition in eq. (Q), it is 
worth noting that only the EM algorithm directly 
optimizes the log-likelihood in eq. (g). This has ob- 
vious advantages if the goal of finding word classes is 
to improve the perplexity of a language model. The 
EM algorithm also handles probabilistic constraints 
in a natural way, allowing words to belong to more 
than one class if this increases the overall likelihood. 

Our approach differs in important ways from the 
use of hidden Markov models (HMMs) for class- 
based language modeling (Jelinek et al., 1992). 



While HMMs also use hidden variables to represent 
word classes, the dynamics are fundamentally dif- 
ferent. In HMMs, the hidden state at time t + 1 is 
predicted (via the state transition matrix) from the 
hidden state at time t. On the other hand, in aggre- 
gate Markov models, the hidden state at time t + 1 
is predicted (via the matrix P(ct+i\u>t)) from the 
word at time t. The state-to-state versus word-to- 
state dynamics lead to different learning algorithms. 
For example, the Baum- Welch algorithm for HMMs 
requires forward and backward passes through each 
training sentence, while the EM algorithm we use 
does not. 

We trained aggregate Markov models with 2, 4, 
8, 16, and 32 classes. Figure [l] shows typical plots 
of the training and test set perplexities versus the 
number of iterations of the EM algorithm. Clearly, 
the two curves are very close, and the monotonic 
decrease in test set perplexity strongly suggests lit- 
tle if any overfitting, at least when the number of 
classes is small compared to the number of words in 
the vocabulary. Table |l] shows the final perplexities 
(after thirty-two iterations of EM) for various ag- 
gregate Markov models. These results confirm that 
aggregate Markov models are intermediate in accu- 
racy between unigram (C = 1) and bigram (C = V) 
models. 

The aggregate Markov models were also observed 
to discover meaningful word classes. Table || shows, 
for the aggregate model with C = 32 classes, the 
most probable class assignments of the three hun- 
dred most commonly occurring words. To be precise, 
for each class c* , we have listed the words for which 
c* = argmax c P(c\w). Figure || shows a histogram of 
the winning assignment probabilities, max c P(c\w), 
for these words. Note that the winning assignment 
probabilities are distributed broadly over the inter- 
val [^j,1]- This demonstrates the utility of allowing 
"soft" membership classes: for most words, the max- 
imum likelihood estimates of P(c\w) do not corre- 
spond to a winner-take-all assignment, and therefore 
any method that assigns each word to a single class 
("hard" clustering), such as those used by Brown et 
al. (1992) or Ney, Essen, and Kneser (1994), would 
lose information. 

We conclude this section with some final com- 
ments on overfitting. Our models were trained by 
thirty-two iterations of EM, allowing for nearly com- 
plete convergence in the log-likelihood. Moreover, 
we did not implement any flooring constraints^ on 

1 lt is worth noting, in this regard, that individual 
zeros in the matrices P(iU2\c) and P(c\wi) do not nec- 
essarily give rise to zeros in the matrix P(u>2\wi), as 
computed from eq. (fil). 




Figure 1: Plots of (a) training and (b) test perplexity versus number of iterations of the EM algorithm, for 
the aggregate Markov model with C = 32 classes. 
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ago day earlier Friday Monday month quarter 
reported said Thursday trading Tuesday 
Wednesday (...) 


20 


did {") c> 


21 


but called San (:) (start-of-sentence) 


22 


bank board chairman end group members 
number office out part percent price prices rate 
sales shares use 


3 


even get to 


4 


based days down home months up work years 

{%> 


23 


a an another any dollar each first good her his its 
my old our their this 


5 


those (,) ( — ) 


6 


(•> <?) 


24 


long Mr. year 


7 


eighty fifty forty ninety seventy sixty thirty 
twenty {() (■) 


25 


business California case companies corporation 
dollars incorporated industry law money 
thousand time today war week {)} (unknown) 


8 


can could may should to will would 


9 


about at just only or than (&} {;) 


26 


also government he it market she that there 
which who 


10 


economic high interest much no such tax united 
well 


27 


A. B. C. D. E. F. G. I. L. M. N. P. R. S. T. U. 


28 


both foreign international major many new oil 
other some Soviet stock these west world 


11 


president 


12 


because do how if most say so then think very 
what when where 


29 


after all among and before between by during for 
from in including into like of off on over since 
through told under until while with 


13 


according back expected going him plan used way 


15 


don't I people they we you 


30 


eight fifteen five four half last next nine oh one 
second seven several six ten third three twelve 
two zero {-) 


16 


Bush company court department more officials 
police retort spokesman 


17 


former the 


31 


are be been being had has have is it's not still 
was were 


18 


American big city federal general house military 
national party political state union York 


32 


chief exchange news public service trade 



Table 2: Most probable assignments for the 300 most frequent words in an aggregate Markov model with 
C = 32 classes. Class 14 is absent because it is not the most probable class for any of the selected words.) 



the probabilities P(c\wi) or P(w2\c). Nevertheless, 
in all our experiments, the ML aggregate Markov 
models assigned non-zero probability to all the bi- 
grams in the test set. This suggests that for large 
vocabularies there is a useful regime 1 <C C -C V 
in which aggregate models do not suffer much from 
ovcrfitting. In this regime, aggregate models can be 
relied upon to compute the probabilities of unseen 
word combinations. We will return to this point in 



Section 4, when we consider how to smooth n-gram 
language models. 

3 Mixed-order Markov models 

One of the drawbacks of n-gram models is that their 
size grows rapidly with their order. In this section, 
we consider how to make predictions based on a con- 
vex combination of pairwise correlations. This leads 



to language models whose size grows linearly in the 
number of words used for each prediction. 

For each k > 0, the skip-k transition matrix 
M(wt-hi w t) predicts the current word from the 
fcth previous word in the sentence. A mixed-order 
Markov model combines the information in these 
matrices for different values of k. Let m denote 
the number of bigram models being combined. The 
probability distribution for these models has the 
form: 



P(w t \w t - 



,W t - m ) 



(6) 



k-1 



Xk{w t -k) M k (w t -k,w t ) 17 [1 - ^j( w t-j) 

3 = 1 



fe=l 



The terms in this equation have a simple interpreta- 
tion. The V x V matrices Mk(w,w') in eq. (^) de- 
fine the skip-fc stochastic dependency of w' at some 
position t on w at position t — k; the parameters 
X k (w) are mixing coefficients that weight the predic- 
tions from these different dependencies. The value of 
X k (w) can be interpreted as the probability that the 
model, upon seeing the word Wt- k , looks no further 
back to make its prediction (Singer, 1996). Thus the 
model predicts from Wt—i with probability Xi(wt-i), 
from wt-2 with probability [1 — Ai(uJt-i)]A2(ti>t-2)) 
and so on. Though included in eq. (|J) for cosmetic 
reasons, the parameters X m (w) are actually fixed to 
unity so that the model never looks further than m 
words back. 

We can view eq. @ as a hidden variable model. 
Imagine that we adopt the following strategy to pre- 
dict the word at time t. Starting with the previous 
word, we toss a coin (with bias Xi(wt-i)) to see if 
this word has high predictive value. If the answer 
is yes, then we predict from the skip-1 transition 
matrix, Mi(w f _i,w f ). Otherwise, we shift our at- 
tention one word to the left and repeat the process. 
If after m — 1 tosses we have not settled on a pre- 
diction, then as a last resort, we make a prediction 
using M m {wt- m ,Wt). The hidden variables in this 
process are the outcomes of the coin tosses, which 
are unknown for each word Wt-k- 

Viewing the model in this way, we can derive an 
EM algorithm to learn the mixing coefficients X k (w) 
and the transition matrices^ M k (w, w'). The E-step 
of the algorithm is to compute, for each word in the 



2 Note that the ML estimates of M k (w, w') do not 
depend only on the raw counts of fc-separated bigrams; 
they are also coupled to the values of the mixing coef- 
ficients, Xk(w). In particular, the EM algorithm adapts 
the matrix elements to the weighting of word combina- 
tions in eq. (^]). The raw counts of fc-separated bigrams, 
however, do give good initial estimates. 




iteration of EM 



Figure 3: Plot of (training set) perplexity versus 
number of iterations of the EM algorithm. The re- 
sults are for the m — 4 mixed-order Markov model. 

training set, the posterior probability that it was 
generated by M k (w t - k ,w t )- Denoting these poste- 
rior probabilities by (j>k(t), we have: 



(7) 



Xk(w t 



-k)M k (w t ^ k ,w t )Yl k j= l[l- 



•Aj(M?t-j)] 



P(w t \wt-l,Wt-2i ■ ■ • , Wt-m) 

where the denominator is given by eq. (||). The 
M-step of the algorithm is to update the parame- 
ters Xk{w) and M k (w,w') to reflect the statistics in 
eq. (0) . The updates for mixed-order Markov models 
are given by: 



Xk(w) 
M k (w,w') 



J2 t d{ w , w t-k)<fik (t) 

J2t S ( w > w t -k)S(w', w t )4> k {t) 

J2 t d( w > w t-k)4>k{t) 



(8) 
, (9) 



where the sums are over all the sentences in the 
training set, and 5(w, w') — 1 iff w — w' . 

We trained mixed-order Markov models with 2 < 
m < 4. Figure || shows a typical plot of the train- 
ing set perplexity as a function of the number of 
iterations of the EM algorithm. Table |j| shows the 
final perplexities on the training set (after four iter- 
ations of EM). Mixed-order models cannot be used 
directly on the test set because they predict zero 
probability for unseen word combinations. Unlike 
standard n-gram models, however, the number of 
unseen word combinations actually decreases with 
the order of the model. The reason for this is that 
mixed-order models assign finite probability to all n- 
grams wiw 2 ■ ■ ■ w n for which any of the ^-separated 
bigrams w k w n are observed in the training set. To 
illustrate this point, Table || shows the fraction of 
words in the test set that were assigned zero proba- 
bility by the mixed-order model. As expected, this 
fraction decreases monotonically with the number of 
bigrams that are mixed into each prediction. 
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train 


missing 


1 


123.2 


0.045 


2 


89.4 


0.014 


3 


77.9 


0.0063 


4 


72.4 


0.0037 



Table 3: Results for ML mixed-order models; m de- 
notes the number of bigrams that were mixed into 
each prediction. The first column shows the per- 
plexities on the training set. The second shows the 
fraction of words in the test set that were assigned 
zero probability. The case to — 1 corresponds to a 
ML bigram model. 

Clearly, the success of mixed-order models de- 
pends on the ability to gauge the predictive value 
of each word, relative to earlier words in the same 
sentence. Let us see how this plays out for the 
second-order (to = 2) model in Table ||. In this 
model, a small value for Xi(w) indicates that the 
word w typically carries less information that the 
word that precedes it. On the other hand, a large 
value for \i(w) indicates that the word w is highly 
predictive. The ability to learn these relationships 
is confirmed by the results in Table |i[ Of the three- 
hundred most common words, Table ^ shows the 
fifty with the lowest and highest values of \i(w). 
Note how low values of Xi(w) are associated with 
prepositions, mid-sentence punctuation marks, and 
conjunctions, while high values are associated with 
"contentful" words and end-of-sentence markers. (A 
particularly interesting dichotomy arises for the two 
forms "a" and "an" of the indefinite article; the lat- 
ter, because it always precedes a word that begins 
with a vowel, is inherently more predictive.) These 
results underscore the importance of allowing the 
coefficients Ai(uj) to depend on the context w, as 
opposed to being context-independent (Ney, Essen, 
and Kneser, 1994). 

4 Smoothing 

Smoothing plays an essential role in language models 
where ML predictions are unreliable for rare events. 
In n-gram modeling, it is common to adopt a re- 
cursive strategy, smoothing bigrams by unigrams, 
trigrams by bigrams, and so on. Here we adopt a 
similar strategy, using the (to — l)th mixed-order 
model to smooth the mth one. At the "root" of 
our smoothing procedure, however, lies not a uni- 
gram model, but an aggregate Markov model with 
C > 1 classes. As shown in Section ||, these models 
assign finite probability to all word combinations, 



O.K AiH < 0.7 
(-) and of <") or <;) to {,) (&> by with S. from 
nine were for that eight low seven the (() (:) six 
are not against was four between a their two 
three its {unknown} B. on as is { — ) five ()} into 
CM. her him over than A. 

0.96 < Xi{w) < 1 
officials prices which go way he last they earlier 
an Tuesday there foreign quarter she former 
federal don't days Friday next Wednesday {%) 
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Tabic 4: Words with low and high values of \i(w) 
in an to = 2 mixed order model. 
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validation 


test 


unseen 


1 


163.615 


167.112 


293175 


2 


162.982 


166.193 


259360 


4 


161.513 


164.363 


200067 


8 


161.327 


164.104 


190178 


16 


160.034 


162.686 


164673 


32 


159.247 


161.683 


150958 



Table 5: Perplexities of bigram models smoothed by 
aggregate Markov models with different numbers of 
classes (C). 

even those that are not observed in the training set. 
Hence, they can legitimately replace unigrams as the 
base model in the smoothing procedure. 

Let us first examine the impact of replacing uni- 
gram models by aggregate models at the root of the 
smoothing procedure. To this end, a held-out inter- 
polation algorithm (Jelinek and Mercer, 1980) was 
used to smooth an ML bigram model with the aggre- 
gate Markov models from Section 0. The smoothing 
parameters, one for each row of the bigram transi- 
tion matrix, were estimated from a validation set the 
same size as the test set. Table || gives the final per- 
plexities on the validation set, the test set, and the 
unseen bigrams in the test set. Note that smooth- 
ing with the C = 32 aggregate Markov model has 
nearly halved the perplexity of unseen bigrams, as 
compared to smoothing with the unigram model. 

Let us now examine the recursive use of mixed- 
order models to obtain smoothed probability esti- 
mates. Again, a held-out interpolation algorithm 
was used to smooth the mixed-order Markov models 
from Section ||. The mth mixed-order model had 
mV smoothing parameters a^iw), corresponding to 
the V rows in each skip-fc transition matrix. The 



to 


validation 


test 


1 


160.1 


161.3 


2 


135.3 


136.9 


3 


131.4 


133.5 


4 


131.2 


133.7 



backoff 


test 


unseen 


baseline 
mixed 


95.2 
79.8 


2799 
1363 



Table 6: Perplexities of smoothed mixed-order mod- 
els on the validation and test sets. 



TOth mixed-order model was smoothed by discount- 
ing the weight of each skip-fc prediction, then fill- 
ing in the leftover probability mass by a lower-order 
model. In particular, the discounted weight of the 
skip-fc prediction was given by 

fc-i 

[1 - a k {w t -k)}h{m-k) - Xjiwt-j)] , (10) 



leaving a total mass of 



fe=i 



fc-i 

a k (w t ~k)Xk(w t -k) Hi 1 ~ Xj(w t -j)] (11) 
i=i 



for the (m — l)th mixed-order model. (Note that 
the to = 1 mixed-order model corresponds to a ML 
bigram model.) 

Table || shows the perplexities of the smoothed 
mixed-order models on the validation and test sets. 
An aggregate Markov model with C = 32 classes 
was used as the base model in the smoothing proce- 
dure. The first row corresponds to a bigram model 
smoothed by a aggregate Markov model; the second 
row corresponds to an to — 2 mixed-order model, 
smoothed by a ML bigram model, smoothed by an 
aggregate Markov model; the third row corresponds 
to an m — 3 mixed-order model, smoothed by a 
m = 2 mixed-order model, smoothed by a ML bi- 
gram model, etc. A significant decrease in perplex- 
ity occurs in moving to the smoothed m = 2 mixed- 
order model. On the other hand, the difference in 
perplexity for higher values of to is not very dra- 
matic. 

Our last experiment looked at the smoothing of 
a trigram model. Our baseline was a ML trigram 
model that backed off to bigrams (and when nec- 
essary, unigrams) using the Katz backoff procedure 
(Katz, 1987). In this procedure, the predictions of 
the ML trigram model are discounted by an amount 
determined by the Good- Turing coefficients; the left- 
over probability mass is then filled in by the backoff 

3 We used a backoff procedure (instead of interpo- 
lation) to avoid the estimation of trigram smoothing 
parameters. 



Table 7: Perplexities of two smoothed trigram mod- 
els on the test set and the subset of unseen word 
combinations. The baseline model backed off to bi- 
grams and unigrams; the other backed off to the 
to = 2 model in Table pi 



model. We compared this to a trigram model that 
backed off to the to = 2 model in Table ||. This was 
handled by a slight variant of the Katz procedure 
(Dagan, Pereira, and Lee, 1994) in which the mixed- 
order model substituted for the backoff model. 

One advantage of this smoothing procedure is that 
it is straightforward to assess the performance of dif- 
ferent backoff models. Because the backoff models 
are only consulted for unseen word combinations, 
the perplexity on these word combinations serves as 
a reasonable figure-of-merit. 

Table ^ shows those perplexities for the two 
smoothed trigram models (baseline and backoff). 
The mixed-order smoothing was found to reduce 
the perplexity of unseen word combinations by 51%. 
Also shown in the table are the perplexities on the 
entire test set. The overall perplexity decreased 
by 16% — a significant amount considering that only 
24% of the predictions involved unseen word com- 
binations and required backing off from the trigram 
model. 

The models in Table |^ were constructed from all 
n-grams (1 < n < 3) observed in the training data. 
Because many n-grams occur very infrequently, a 
natural question is whether truncated models, which 
omit low- frequency n-grams from the training set, 
can perform as well as untruncated ones. The ad- 
vantage of truncated models is that they do not need 
to store nearly as many non-zero parameters as un- 
truncated models. The results in Table || were ob- 
tained by dropping trigrams that occurred less than 
t times in the training corpus. The t = 1 row cor- 
responds to the models in Table [?]. The most in- 
teresting observation from the table is that omitting 
very low-frequency trigrams does not decrease the 
quality of the mixed-order model, and may in fact 
slightly improve it. This contrasts with the standard 
backoff model, in which truncation causes significant 
increases in perplexity. 



t 


baseline 


mixed 


trigrams(x 


10 b ) 


missing 


1 


95.2 


79.8 


25.4 




0.24 


2 


98.6 


78.3 


6.1 




0.32 


3 


101.7 


79.6 


3.3 




0.36 


4 


104.2 


81.1 


2.3 




0.38 


5 


106.2 


82.4 


1.7 




0.41 



Table 8: Effect of truncating trigrams that occur 
less than t times. The table shows the baseline and 
mixed-order perplexities on the test set, the num- 
ber of distinct trigrams with t or more counts, and 
the fraction of trigrams in the test set that required 
backing off. 

5 Discussion 

Our results demonstrate the utility of language mod- 
els that are intermediate in size and accuracy be- 
tween different order n-gram models. The two 
models considered in this paper were hidden vari- 
able Markov models trained by EM algorithms for 
maximum likelihood estimation. Combinations of 
intermediate-order models were also investigated by 
Rosenfeld (1996). His experiments used the 20,000- 
word vocabulary Wall Street Journal corpus, a pre- 
decessor of the NAB corpus. He trained a maximum- 
entropy model consisting of unigrams, bigrams, tri- 
grams, skip-2 bigrams and trigrams; after selecting 
long-distance bigrams (word triggers) on 38 million 
words, the model was tested on a held-out 325 thou- 
sand word sample. Rosenfeld reported a test-set 
perplexity of 86, a 19% reduction from the 105 per- 
plexity of a baseline trigram backoff model. In our 
experiments, the perplexity gain of the mixed-order 
model ranged from 16% to 22%, depending on the 
amount of truncation in the trigram model. 

While Rosenfeld's results and ours are not di- 
rectly comparable, both demonstrate the utility of 
mixed-order models. It is worth discussing, how- 
ever, the different approaches to combining infor- 
mation from non-adjacent words. Unlike the max- 
imum entropy approach, which allows one to com- 
bine many non-independent features, ours calls for 
a careful Markovian decomposition. Rosenfeld ar- 
gues at length against naive linear combinations in 
favor of maximum entropy methods. His arguments 
do not apply to our work for several reasons. First, 
we use a large number of context-dependent mixing 
parameters to optimize the overall likelihood of the 
combined model. Thus, the weighting in eq. (](]) en- 
sures that the skip-/c predictions are only invoked 
when the context is appropriate. Second, we adjust 
the predictions of the skip-/c transition matrices (by 



EM) so that they match the contexts in which they 
are invoked. Hence, the count-based models are in- 
terpolated in a way that is "consistent" with their 
eventual use. 

Training efficiency is another issue in evaluating 
language models. The maximum entropy method 
requires very long training times: e.g., 200 CPU- 
days in Rosenfeld's experiments. Our methods re- 
quire significantly less; for example, we trained the 
smoothed m — 2 mixed-order model, from start to 
finish, in less than 12 CPU- hours (while using a 
larger training corpus). Even accounting for differ- 
ences in processor speed, this amounts to a signifi- 
cant mismatch in overall training time. 

In conclusion, let us mention some open problems 
for further research. Aggregate Markov models can 
be viewed as approximating the full bigram tran- 
sition matrix by a matrix of lower rank. (From 
eq. (Q) , it should be clear that the rank of the class- 
based transition matrix is bounded by the num- 
ber of classes, C.) As such, there are interesting 
parallels between Expectation-Maximization (EM), 
which minimizes the approximation error as mea- 
sured by the KL divergence, and singular value de- 
composition (SVD), which minimizes the approxi- 
mation error as measured by the Li norm (Press 
et al, 1988; Schutze, 1992). Whereas SVD finds a 
global minimum in its error measure, however, EM 
only finds a local one. It would clearly be desirable 
to improve our understanding of this fundamental 
problem. 

In this paper we have focused on bigram models, 
but the ideas and algorithms generalize in a straight- 
forward way to higher-order n-grams. Aggregate 
models based on higher-order n-grams (Brown et al., 
1992) might be able to capture multi-word struc- 
tures such as noun phrases. Likewise, trigram-based 
mixed-order models would be useful complements to 
4-gram and 5-gram models, which are not uncom- 
mon in large-vocabulary language modeling. 

A final issue that needs to be addressed is 
scaling — that is, how the performance of these mod- 
els depends on the vocabulary size and amount 
of training data. Generally, one expects that the 
sparser the data, the more helpful are models that 
can intervene between different order n-grams. Nev- 
ertheless, it would be interesting to see exactly how 
this relationship plays out for aggregate and mixed- 
order Markov models. 
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